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Abstract. Non-steady and eruptive phenomena in quasars are thought to be associ- 
ated with the Eddington or super-Eddington luminous stage. Although there is no lack 
in hypotheses about the total duration of such a stage, the latter remains essentially 
unknown. We calculate the duration of quasar luminous phase in dependence upon 
the initial mass of a newborn massive black hole (MBH) by comparing the observed 
luminosity- and redshift distributions of quasars with mass distribution of the central 
MBHs in normal galactic nuclei. It is assumed that, at the quasar stage, each MBH 
goes through a single (or recurrent) phase(s) of accretion with, or close to, the Ed- 
dington luminosity. The mass distributions of quasars is found to be connected with 
that of MBHs residing in normal galaxies by a one-to-one corrrespondence through the 
entire mass range of the inferred MBHs if the accretion efficiency of mass-to-energy 
transformation r\ ~ 0.1. 



INTRODUCTION 

An approximate relationship Mh — (0.003 — .006)M& between the central MBH 
mass, Mh, and that of the galactic bulge, M&, has been established for a few dozen 
of galaxies, both nearby and more distant ones [1], [2]. A relationship between 
absolute magnitudes of quasars and their host galaxies found in [3] is reduced to 
the MBH to bulge mass relation in galaxies provided that [4]: (i) the central MBH 
shines at or near to the Eddington luminosity and (ii) the host galaxy undergoes 
through a starburst episode. This correlation, coupled with the known luminosity 
function of galaxies, can serve to obtain [5] the MBH mass distribution 1 (M/ l )c?M/ l . 
The history of matter accretion onto a central MBH thought to serve as a source 
of quasar activity is linked to the present observable properties of each individual 
quasar, such as its luminosity, variability, and emission spectrum. If the bolometric 
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luminosity of a quasar comprises a fraction A of the Eddington luminosity, i.e. 
A = L/Le, Le = 4:iiGMhm p c/o-T, the underlying accretion is accompanied by 
an exponential growth of the MBH mass with the characteristic time t E = 4.5 • 
10 8 ?7/A yrs, where 77 is the accretion efficiency of mass-to-energy transformation. 
The crucial problem is the duration, t q , of such a nearly Eddington accretion phase. 
Usually an effective t q , the same for the entire black hole mass range M h ~ 10 6 — 
1O 1O M , is calculated by comparing the global number density of normal galaxies 
and quasars and is found to be t q = 10 6 — 5 • 10 8 years in [6]. Meanwhile the recent 
data on mass distribution of MBHs in galaxies provide an opportunity to solve this 
problem in a more detailed way, viz., to calculate the dependence of t q upon M h , 
which is the major aim of this paper. We shall explore whether the distribution 
functions of quasars and MBHs in normal galaxies are consistent with each other, 
and we will do this locally in the vicinity of each mass. 



THE EDDINGTON LUMINOSITY PHASE 

It would be reasonable to assume that the duration of the Eddington phase t q 
depends on the initial mass of a newborn MBH or, in other words, on the initial 
luminosity, Lj, of the quasar: t q = t q (Li). For simplicity, the transition to and out 
of the Eddington phase is supposed to occur instantaneously: 

0, if t < U; 

Li€sxp[(t-ti)/t E ], if U < t < U + t q (Li); (1) 
0, if t >t t + t q {L t ), 

where U is the instant of the MBH formation. Along with the distribution function 
of MBHs in the galactic nuclei, 0i, we use the observed distribution of quasars in 
absolute magnitude Mb and redshift z < z e ~ 3, ^(Mb, z)dMsdz, taken from [7]. 
The balance equation is given by 

00 . . x 

— Jdzil + zy^MMBiL)^)^*-^- J MM h )dM h , (2) 

M(L) 

where X = M(L)exp[t q (L)/tE] with M(L) determined by equation L = \L E and 
a relationship t = t /(l + z) 3 ^ 2 for the flat cosmological model is used. While 
obtaining Eq. (2), we have also taken into account that t E very slowly varies with 
respect to <f>i(Mh). Eq. (2) defines t q as an implicit function of L iy which can 
be translated into a relationship between t q and the initial BH mass Mj by using 
equation L = \L E . 



RESULTS 



Numerical solution of Eq. (2) is found by adopting the mass distribution of 
MBHs (f>i(M h )dM h from [5], derived with the use of three relationships, viz., (i) 
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FIGURE 1. The ratio t q (Mi)/t E as a function of the initial MBH mass M t . Curves labeled 1 
to 6 are based on distribution A for rj =10 _1 , 1CP 2 , 1CP 3 , 1CP 4 , 10~ 5 , and 1CP 6 , respectively. 
Curves 7 and 8 are based on distribution B for rj =10^ 1 and 10 -2 , respectively. 

a correlation \og(M h ) = log(M 6 ) — 2.6 ± 0.3 between the MBH mass M h and 
the bulge mass M b ; (ii) the mass-luminosity relation for galaxies, and (iii) the 
Schechter luminosity function. In Fig. 1 and Fig. 2, we employ two somewhat 
different distributions in MBH mass from [5] and name them 'distribution A' and 
'distribution B', which correspond to the power-law and log-Gaussian shape of 
dispersion, respectively. Fig. 1 presents the results of our numerical computation 
of the ratio t q (Mj)/tE for different values of 77 and MBH mass distributions A and 
B. It should be noted that if rj < 0.1, the solution exists not for all values of Mj. 
The domain where the solution exists is determined from the condition that the 
r.h.s of Eq. (2) exceeds its l.h.s. if one puts X = +00. For those Mj which lead 
to an opposite condition, the number of galactic nuclei with MBHs is not enough 
to explain, in the framework of our model, the distribution function of quasars in 
Mb and z, even if these MBHs stay in an active quasar state during the maximum 
possible time t q ~ 3t^. The solution only exists at r\ > 7 ■ 10~ 7 for BH mass 
distribution A and at r\ > 6 • 10~ 3 for distribution B. A single-valued mapping 
Mi — > Mh breaks up on the left end of curves 2 to 6. 

Fig. 1 demonstrates the main result of this work: distributions of MBHs and 
quasars in mass are connected by one-to-one correspondence through the whole 
range of the observed masses only for r\ ~ 0.1, both for the distribution A (curve 
1) and B (curve 7). This concordance breaks down for a certain range of MBH 
masses, viz. the solution of Eq. (2) does not exist if i] <C 0.1. Nevertheless, the 
jumps of the r\ value are not excluded on the boundary of the domain, where the 



10 6 10 7 10 8 10 9 10 1 ° 10 11 10 12 

M i / M 

FIGURE 2. The accreted mass, Mf — Mi, as a function of the initial mass Mi. The curves are 
labeled in the same way as in Fig. 1. 

solution exists. Similar jumps seem to be quite natural if MBHs in the different 
mass ranges are formed by different ways (e.g., by collapse of massive gas clouds, 
stellar clusters, etc., see review by Rees 1984) and so there are various accretion 
regimes with different values of rj. If such jumps indeed take place, transitions 
between the curves of each of distributions A and B are possible. These transitions 
must be smoothed because the MBHs formed by different ways would coexist in 
some mass range(s). The most probable value of r] established above is r) ~ 0.1, 
which corresponds to curves 1 and 7 in Fig. 1. For both these curves, the relation- 
ship t q < t E is carried out and therefore the BH mass growth is not substantial 
- it generally does not exceed a value comparable to the initial BH mass, and for 
M h > 10 7 M it is negligible. 
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